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We provide a general framework to describe cooling of a micromechanical oscillator to its quantum 
ground state by means of radiation-pressure coupling with a driven optical cavity. We apply it to 
two experimentally realized schemes, back-action cooling via a detuned cavity and cold-damping 
quantum-feedback cooling, and we determine the ultimate quantum limits of both schemes for the 
full parameter range of a stable cavity. While both allow to reach the oscillator's quantum ground 
state, we find that back-action cooling is more efficient in the good cavity limit, i.e. when the cavity 
bandwidth is smaller than the mechanical frequency, while cold damping is more suitable for the bad 
cavity limit. The results of previous treatments are recovered as limiting cases of specific parameter 
regimes. 
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INTRODUCTION 

Cooling of mechanical resonators close to their quan- 
tum ground state has become an important topic for vari- 
ous fields of physics, such as ultra-high precision measure- 
ments in , the detection of gravitational waves 0] , and the 
study of the transition between classical and quantum 
behavior of a mechanical system [^. It is also a prereq- 
uisite for any possible use of optomechanical systems for 
quantum information processing 0, Q . Recently, various 
experiments have demonstrated significant cooling of the 
vibrational mode of a mechanical resonator coupled to 
an optical cavity @, S B i, 113, Ei El, [11, H, 111, 13. 
In these experiments cooling has been achieved by ex- 
ploiting in two different ways the radiation-pressure inter- 
action between a mechanical mode and the intra-cavity 
field: i) by back-action, or self-cooling 17[, in which the 
off-resonant operation of the cavity results in a retarded 
back action on the mechanical system and hence in a 
"self" -modification of its dynamics [ol . [Til , [isl . [lil [ist : ii) 
by cold-damping quantum feedback [l9|, |20|, l2ll |. where 
the oscillator position is measured through a phase- 
sensitive detection of the cavity output and the result- 
ing photocurrent is used for a real-time correction of the 
dynamics [1, [H, [H, [H . 

We generalize and extend the previous treatments of 
these schemes [11, [H, [H, [2l|, [H, [Ig to the full parameter 
range of a stable cavity by deriving the quantum steady 
state of the micromechanical oscillator in a linearized 
quantum Langevin equation (QLE) approach. Compar- 
ing the two schemes we find that back-action cooling is 
more efficient in the good cavity limit, i.e. when the cav- 
ity bandwidth is smaller than the mechanical frequency, 
while cold damping is more suitable in the opposite limit 
of a bad cavity. We also show that, contrary to common 
belief, the feedback gain in cold-damping schemes is nec- 



essarily bounded by an upper limit to achieve quantum 
ground-state cooling. 

The paper is organized as follows. In Section II we de- 
scribe the dynamics of the system in terms of linearized 
quantum Langevin equations. In Sec. Ill and IV we eval- 
uate the steady-state energy of the mechanical oscillator 
for the two cases of back-action cooling with a detuned 
cavity and cold-damping feedback cooling. In Section 
V we conclude by comparing in detail the two cooling 
schemes. 



QUANTUM LANGEVIN EQUATIONS FOR THE 
SYSTEM 



We consider a driven optical cavity coupled by radia- 
tion pressure to a micromechanical oscillator. The typ- 
ical experimental configuration is a Fabry-Perot cavity 
with one mirror much lighter than the other (see e.g. 



ll|, [ij, ll^l), but our treatment applies to other 



configurations such as the silica toroidal microcavity of 
Refs. [l3l.[2^. Radiation pressure typically excites several 
mechanical degrees of freedom of the system with differ- 
ent resonant frequencies. However, a single mechanical 
mode can be considered when a bandpass filter in the 
detection scheme is used [11] and coupling between the 
different vibrational modes can be neglected. The Hamil- 
tonian of the system reads 24 1 
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The first term describes the energy of the cavity mode, 



with lowering operator a ([o 



I), cavity frequency 



uic and decay rate k. The second term gives the en- 
ergy of the mechanical mode, modeled as harmonic os- 
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cillator at frequency and described by dimensionless 
position and momentum operators q and p {[q,p] = i). 
The third term is the radiation-pressure coupling of rate 
Go = (ijJcl L)^Jh/mujjn, where m is the effective mass 
of the mechanical mode [2^ , and L is an effective length 
that depends upon the cavity geometry: it coincides with 
the cavity length in the Fabry-Perot case, and with the 
toroid radius in the case of Refs. [3, 22|. The last term 



describes the input driving by a laser with frequency 
Wo I where E is related to the input laser power P by 
\E\ = y/2PK/hjjJa. One can adopt the single cavity mode 
description of Eq. ^ as long as one drives only one cavity 
mode and the mechanical frequency is much smaller 
than the cavity free spectral range FSR ~ c/L. In this 
case, scattering of photons from the driven mode into 
other cavity modes is negligible [2^. 

The dynamics are also determined by the fluctuation- 
dissipation processes affecting both the optical and the 
mechanical mode. They can be taken into account in a 
fully consistent way [2J| by considering the following set 
of nonlinear QLE, written in the interaction picture with 
respect to hhjQO^ a 



q = UJmP, 

p = ~uj„iq - jmP + Goa'^a + ^, 
d — —{K + iAQ)a + iGQaq + E 



(2a) 
(2b) 

2Ka"\ (2c) 



where Ao = lUc — loq. The mechanical mode is affected by 
a viscous force with damping rate 7^ and by a Brownian 
stochastic force with zero mean value ^, that obeys the 
correlation function 24l. |27| 



mm) 



duj 
2^' 



-iu{t — t') 



coth 



2kBT 



(3) 

where fc^ is the Boltzmann constant and T is the temper- 
ature of the reservoir of the micromechanical oscillator. 
The Brownian noise ^(t) is a Gaussian quantum stochas- 
tic process and its non-Markovian nature (neither its cor- 
relation function nor its commutator are proportional to 
a Dirac delta) guarantees that the QLE of Eqs. ([2|) pre- 
serve the correct commutation relations between oper- 
ators during the time evolution 2J|. The cavity mode 
amplitude instead decays at the rate k and is affected by 
the vacuum radiation input noise a"^{t), whose correla- 
tion functions are given by 



and 



(a™(i)a"''^(i')> - [N{i^c) + 1] Sit - t'), (4) 



where N{oJc) — (explfi^c/fesT} — 1) ^ is the equilib- 
rium mean thermal photon number. At optical frequen- 
cies HcUc/kBT ^ 1 and therefore N{u!c) ~ 0, so that only 
the correlation function of Eq. ^ is relevant. 



Cooling of the mechanical oscillator by radiation pres- 
sure can be described in thermodynamical terms in the 
following way. Radiation pressure couples the oscillator 
to the optical cavity mode, which behaves as an effective 
additional reservoir for the oscillator when the cavity is 
appropriately detuned. As a consequence, the effective 
temperature of the mechanical mode will be intermedi- 
ate between the initial reservoir temperature and that of 
the effective optical reservoir, which is in practice equal 
to zero due to the condition N{oJc) — 0. Therefore one 
approaches the mechanical ground state when the cou- 
pling rate to the optical reservoir is much larger than the 
damping rate which gives the coupling to the initial 
reservoir. This explains why significant cooling is ob- 
tained when radiation pressure coupling is strong. It is 
realized when the coupling Go is large, but is more easily 
achieved when the intracavity field is very intense, i.e., for 
high-finesse cavities and enough driving power. In this 
limit (and if the system is stable) the system is character- 
ized by a semiclassical steady state with the cavity mode 
in a coherent state with amplitude as {\cts\ ^ 1), and 
a new equilibrium position for the oscillator, displaced 
by qs- The parameters as and qs are the solutions of 
the nonlinear algebraic equations obtained by factorizing 
Eqs. ([2|) and setting the time derivatives to zero. They 
are given by 



Go\as\' 
E 



(6) 
(7) 



where the latter equation is in fact the nonlinear equation 
determining as, since the effective cavity detuning A, 
including radiation pressure effects, is given by 



A = Ao 



Gl\as 



(8) 



Rewriting each Heisenberg operator of Eqs. ([2]) as the c- 
number steady state value plus an additional fluctuation 
operator with zero mean value, one gets the exact QLE 
for the fluctuations 

6q = LUmSp, 

6p = —uj„iSq — jmSp + Go[asSa'+alSa)+Sa^5a + S,, 
Sa = -{K + iA)Sa + iGoias + Sa)Sq + V2Ka"\ (9) 

We have assumed ja^j ^ 1, therefore one can safely ne- 
glect the nonlinear terms Sa'^Sa and SaSq in the equations 
above and obtains the linearized QLE 

Sq — LOmSp, (10a) 
Sp = -UJrnSq ~ -frnSp + GSX + (10b) 

5X = -k5X + A5Y + V2^X'", (10c) 
dY = -kSY - ASX + GSq + V2^Y'''. (lOd) 
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Here we have chosen the phase reference of the cavity 
field so that as is real and positive, and we have de- 
fined the cavity field quadratures SX = {Sa + Sa^ ) / \/2 
and SY = {Sa — Sa^)/i^/2, and the corresponding Her- 
mitian input noise operators X*" = (a™ + a^"''^)/\/2 
and = (a'" - a'"^^)/iV2. The hnearized QLE show 
that the mechanical mode is coupled to the cavity mode 
quadrature fluctuations by the effective optomechanical 
coupling 



G — Goa^V^ = 



2uJr 



Pk 



(11) 



which can be made very large by increasing the intracav- 
ity amplitude Ug . Notice that together with the condition 
'^m c/L which is required for the single cavity mode 
description, ^ 1 is the only assumption required by 
the present approach. We shall see below that, thanks to 
this fact, the present approach provides a generalization 
of previous treatments of back-action cooling, such as the 
semiclassical treatment of Ref. fll\ and the perturbation 
treatments of 28, 29l. 



DETUNING-INDUCED BACK-ACTION 
COOLING 

We have to evaluate the mean energy of the oscillator 
in the steady state 

fl^n. 



U 



[{Si') + i^p')] 



huJri 



(12) 



and to see if and when it approaches the ground-state 
value hujjn/2. This is equivalent to determine the condi- 
tions under which {Sq'^) ~ {^p'^) — 1/2. The two oscilla- 
tor variances (Sq'^) and {Sp'^) can be obtained by solving 
Eqs. (dni) in the frequency domain and integrating the 
corresponding fluctuation spectrum. One gets 



where the position spectrum is given by 
where 



coth 



2kBT 



is the thermal noise spectrum, 
G^n [A2 



Srp{u),A) 



A)2] [k2 + + A)2] 



(13) 
(14) 

(15) 
(16) 



is the radiation pressure noise spectrum, and 



G^Au 



is the effective susceptibility of the oscillator, modified by 
radiation pressure. The latter can be read as the suscep- 
tibility of an oscillator with effective resonance frequency 
and damping rate given by 



(^) 



G^Auj^{k^-uj^ + A^) 

[k2 + {uj- A)2] [k2 + + A)2] 

26*2 Aw„K 



(18) 



(19) 



[k2 + {uJ~ A)2] [^2 + {u; + A)2j 

The modification of the mechanical frequency due to ra 



coeff/co, 




moo\ A 




FIG. 1: Plot of the effective mechanical frequency of Eq. (llSp 
(a) and of the effective mechanical damping of Eq. p9|l (b) 
versus frequency. Parameter values are tJm/27r = 10 MHz, 

7m/27r — 100 Hz, A = LUm, H = 0.2 UJrn and G — 0.2 LUm 

(dashed line), G = 0.3 com (dotted line) and G = 0.4 ujm 
(full line), which we shall see later correspond to an optimal 
cooling regime. 

diation pressure shown by Eq. ^TE\i is the so-called "op- 
tical spring effect" , which may lead to significant fre- 
quency shifts in the case of low-frequency oscillators, 
such as pendulum modes of suspended mirrors [13] ■ In 
the case of higher resonance frequencies, such as those of 
Refs. [i, [n, 13| where tOm ^ 1 MHz, the optical spring 
term in Eq. (fT5|) does not significantly alter the frequency, 
even for large intracavity power. Here, we shall only 
(17) consider numerical examples with large curm where the 
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frequency is practically unchanged {(^^^^{uj) ~ uJm, see 
Fig-Hh)- In fact, ground state cooling rig// < 1 can be 
approached only if the initial mean thermal excitation 
number n = (expl/kj^/fesT} — 1)""'^ is not prohibitively 
large, and this is possible, even at cryogenic tempera- 
tures, only if LUm is sufficiently large. For positive A and 
for large enough G the effective mechanical damping is in- 
stead significantly increased (see Fig.[T|3). This increase 
is at the basis of the cooling process. In fact, the mechan- 
ical susceptibility at resonance is inversely proportional 
to damping and it is therefore significantly suppressed 
by radiation pressure. As a consequence, the oscillator is 
much less affected by thermal noise and this means cool- 
ing as long as the radiation-pressure noise contribution 
Srp{u!, A) remains small compared to thermal noise term 
Sthi^)i which is verified for not too large G. 

Let us now determine the oscillator mean energy. The 
system reaches a steady state only if it is stable and this is 
satisfied when all the poles of the effective susceptibility 
X'tff lis in the lower complex half-plane. By applying 
the Routh-Hurwitz criterion [sO], we get the following 
two nontrivial stability conditions in the detuned-cavity 



Sl 



27mK I 't^ + (t^m - A)^ + {uJrn + A) 



+7r« [(7m + 2k) (k^ + A^) + 2kuj^ } 
+Alu,^G^ (7„ + 2K,f > 0, 



S2 



G^A > 0. 



(20a) 
(20b) 



When the stability conditions are satisfied, the integrals 
of Eq. (fT3)) for the two variances can be solved exactly. 
However, it is reasonable to simplify the thermal noise 
contribution in Eqs. (fT5| . In fact ksT/H ~ 10^^ s~^ even 
at cryogenic temperatures and therefore is always much 
larger than all the other parameters. At these high values 
of uj the position spectrum is negligible and therefore one 
can safely approximate in the integral 



coth ' 



\2kBT 



2kBT 

huJm 



7™(2n+l). (21) 



Performing the integrals, one gets the final expressions 
for the two variances, which are given by 



{Sq'^) = ^ + bg + dqU, 
i + 6p + dpfi, 



(22) 
(23) 



where 



bp = [si] ^G^nlA^i^-frn + k) + K{"f,nH + n"^ + UJm) - AuJmijm +2k)} 
dp = l-[si]-'^2G'^KUJmA{j,n + 2K), 



bq = [2.si52]"'G'2|2k(A2 + k2) I 



A' + (7™ + k) 

^m (7m + k) i^^m - 2A)} + AG^LUm (7„, + 2k) [A7„ - k (tj,„ - 2A)]} , 
1 + [siS2]-^AG^ [si - 2jmHLjl, {ujI, + 27™^ -I- 4^2) _ Ak'ojI^ (a2 + K^)] 

I 



(24) 
(25) 
(26) 

(27) 



At the ground state both variances are equal to 1 /2 and 
therefore realizing ground state cooling means achieving 
bq, bp, dq, dp 0. Notice that in general {5q^) ^ {Sp^), 
that is, one does not have energy equipartition, as it is 
already shown by the general expression of Eq. p3p . This 
means that in the generic case, the steady state of the sys- 
tem is not, strictly speaking, a thermal equilibrium state 
and this prevents to derive from here an univocally de- 
fined temperature. With this respect, Eq. (fT2|) provides 
a definition of the effective mean excitation number riej f 
from which one can only define an effective temperature 
as Teff = hiUm/ [fcBln(l -|- l/ueff)]- Howcver, in order 
to get to the quantum ground state, both variances have 



to tend to 1/2 and therefore energy equipartition has to 
be satisfied in the optimal regime close to the ground 
state. 

In order to have an intuitive picture and to facilitate 
the com par ison with the recent perturbation treatments 
of Refs. H, 2§|, we consider the expressions of the vari- 
ances in some limiting case of experimental interest. 

It is convenient to introduce the rates [291 



A, 



G^K 



2 + {A±ijj„ 



(28) 



which define the rates at which laser photons are scat- 
tered by the moving oscillator simultaneously with the 
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absorption (Stokes, A+) or emission (anti-Stokes, A-) of 
the oscillator vibrational phonons. For A > one has 
A_ > and a net laser cooling rate 



T^A_ 



A, > 



can be defined, giving the rate at which mechanical en- 
ergy is taken away by the leaking cavity. As a conse- 
quence, the total mechanical damping rate is given by 
7m + r, which is consistent with the expression of the ef- 
fective (frequency-dependent) damping rate of Eq. 
in fact, it is easy to check that 



(29) 



As discussed in the preceding section, ground-state cool- 
ing is achievable when 7,„, the coupling rate with the 
thermal reservoir, is significantly smaller than F, which 
represents the coupling rate of the mechanical oscillator 
with the effective reservoir provided by the damped cav- 
ity mode. For small 7™ the above expressions simplify 
to 



(V) 



(Si') 



1 



7m + F 



1 

7m + F 



A, 



+ 7m" 



F 



(30) 



2k 



where we have defined the coefficients 



b = 



T]A (^2 + A2 
2 (A2 - , 



2 



- oj: 



m 



VA 



= 1- 



+ A2 



(31) 

(32) 

(33) 
(34) 



Note that for positive A, < 77A < 1 due to the stability 
condition of Eq. pIb V Eqs. and (EJ) provide a 

generalization of the results of Ref. |28i, l29|, which are 
reproduced if we take Um ^ "■7m , G and k ^ 7„j , G 
(assumed in Refs. [H il) in Eqs. ([201) and dST])- In 
these limits a, ?7a — > 1, F/k — > and therefore 



(Sp^) 0, (Sq^) = ne// + 1/2, 



where 



7mn -I- A+ 
7m + F 



(35) 



(36) 



is the basic resuh of Ref. [28| (see Eq. (9)) and Ref. |29| 
(see Eq. (5) and its derivation). This result has been also 
obtained in Ref. 31| through an approximate treatment 
of the exact integrals of Eq. (fTS]) . 

The approximate expression of Eq. ([36|) suggests how 
the system parameters can be chosen in order to mini- 
mize the oscillator energy close to the ground-state value. 



In order to reach the condition Uef / < 1 one needs a large 
F/7m while keeping A^/T small. This is done by match- 
ing A — LOm, which optimizes the energy transfer from 
the mechanical mode to the anti-Stokes sideband. Fur- 
ther optimization requires a large value of G, which is 
however constrained by the stability condition t^a > 0. 
When 77A — + 0, (Sq^) becomes too large, while the ground 
state, where (Sq^) = (Sp^) = 1/2, is approached for 
?7A 1 (see Eq. ([31])). Finally the optimal ratio k/uj„i is 
determined by the fact that the laser-cooling rate F has 
to be large but still smaller than k because the cavity 
response time has to be remain shorter than that 

of the oscillator F~^. We find that, within a set of 
experimentally achievable parameters, the best cooling 
regime {rieff ~ 0.1), is obtained in the good cavity limit 
condition ^/wm — 0.2 (see Fig. [2t , which is close to the 
value of 1 / V32 suggested in [29| . 
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FIG. 2: (Color online) a) Effective mean vibrational number 
rie// versus A/curn and ft/uim around the optimal ground- 
state cooling regime for ajm/Syr = 10 MHz, 7m/27r = 100 Hz, 
m = 250 ng, a cavity of length L = 0.5 mm driven by a laser 
with power P = 50 mW, and wavelength 1064 nm. The os- 
cillator reservoir temperature is T = 0.6 K, corresponding to 
n ~ 1250. The minimum value rieff — 0.1 corresponds to an 
effective temperature Te// ~ 0.2 mK. (b) Effective mean vi- 
brational number n^ff versus K/um and the normalized power 
P/Po, {Pq ~ 50 mW) at the fixed, optimal value for the de- 
tuning, A — LOrn- The other parameters are the same as in 
(a). 
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GROUND-STATE COOLING WITH COLD 
DAMPING 



so that the estimated phase quadrature SY^^* (s) is given 
by [331] 



An alternative way of cooHng the oscillator by over- 
damping, proposed in [l9| and experimentally realized in 
[6, 10, 12,25.], is to use quantum feedback and specifically 
cold damping 0, 21, HI]- This technique is based on 
the application of a negative derivative feedback, which 
increases the damping of the system without increasing 
the thermal noise. The oscillator position is measured by 
means of a phase-sensitive detection of the cavity out- 
put, which is then fed back to the oscillator by applying 
a force whose intensity is proportional to the time deriva- 
tive of the output signal, and therefore to the oscillator 
velocity. One measures the phase quadrature Y, whose 
Fourier transform, according to Eqs. pO|) . is given by 



SY{uj) 



G{k — itu) 



Sq{uj) 



noise terms. 



(37) 



This shows that the highest sensitivity for position mea- 
surements is achieved for a resonant cavity, A = and 
in the large cavity bandwidth limit k 3> Wm,7m, i.e., 
when the cavity mode adiabatically follows the oscillator 
dynamics 6Y{uj) ~ [G / K)5q{Lo). Therefore the QLE for 
the cold-damping scheme coincide with those of Eqs. (fTT 
with A = 0, but with an additional feedback force 



5q 
6p 



SX 
SY 



-UrnSq - JrnSp + GSX + ^, 

dsg{t- s)6Y^'\s), 



-n5X + V2nX'", 
-k6Y + GSq + V2^y* 



(38a) 
(38b) 



(38c) 
(38d) 



SY'^'^t) 



r°"*(i) 



SY{t) 



y™(i) + y^T]-^ - Wit) 



(41) 

After solving the QLE of Eqs. ([55)1 by Fourier transform, 
one finds that the two oscillator variances in the cold- 
damping case are given again by Eqs. (jl3p . but with a 
different position spectrum S'^fiuS), that can be expressed 



ST{^) = \xl)M\nSth{Lo) + 5.p(c.,0) + Ss^{u)\. (42) 

SthW) is given by Eq. p5)) . S'rp(ci;, 0) is obtained by 
choosing A = in Eq. (fT6|) . and one has the additional 
contribution 



5/b(w) 



(43) 



due to the measurement noise that is fed back into the 
oscillator dynamics by the cold-damping feedback loop 
(^((jj) is the Fourier transform of g{t)). Finally, 



g{uj)Gujr, 



K — lUJ 



n -1 



(44) 



is the effective susceptibility of the oscillator in the cold- 
damping scheme, which depends upon the explicit form 
of the feedback transfer function g{u>). The simplest 
choice, corresponding to a standard derivative high-pass 
filter, is 



9{^) = 



1 — ILU jiO fb ' 



(45) 



which means choosing 



In Eq. ([38b ) g{t) is a causal kernel, proportional to a 
derivative of a Dirac delta in the ideal derivative feedback 
limit, and SY'^''*{s) is the estimated intra-cavity phase 
quadrature that is obtained from the measurement of the 
output quadrature F°"*(t) as follows. The usual input- 
output relation 



SY°^t{t) = V2^SY{t) - r"(i) 



(39) 



can be generalized to the case of a non-unit detection ef- 
ficiency, by modeling a detector with quantum efficiency 
ry with an ideal detector preceded by a beam splitter with 
transmissivity y/rj, hence mixing the incident field with 
an uncorrelated vacuum field Y^ (t) [2^ . The generalized 
input-output relation then reads 



r°"*(t) = v?^ V2^^r(t)-r"(i) - 



(40) 



g{t) = 9cd-j^ [0{t)ujfbe~ 



ulfbt] 



so that 



plays the role of the time delay of the 



feedback loop, and g^d > is the feedback gain. The 
ideal derivative limit is obtained for cofb — *■ cxi, implying 
g{uj) = —iujgcd and therefore g{t) = gcd5'{t). The cold- 
damping susceptibility of Eq. ([44)1 can be read again as 
the susceptibility of an oscillator with effective resonance 
frequency and damping rate given by 



UJ. 



eff.c. 



gcdGujjnUJfbUj'^^K + UJfb) 
(k2 +t^2)(^2^ +^2) 



(46) 



^-ff^cd( , gcdGcj^ujfbjtiUJfb-u;^) 

The frequency dependence of the effective resonance fre- 
quency and damping depends upon the specific form of 
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half-plane we get one nontrivial stability condition 

Scd = [imKUJfb + gcdGuJniLOfb + ^mi^ + ^fb)] 

X [{k + 7,„)(k + UJfb){lrn + i^fb) + Jm(^rn - 9cdGuj.rnOJ fb] 

-KUJ^UJfb{K + Jm+ l^Sbf > 0. (48) 



This stability condition is always satisfied in the ideal 
derivative feedback limit (tOfb — > oo) and in the ideal 
adiabatic limit k oo. However, for finite values of LUfh 
and K, Eq. implies an upper limit to the feedback 
gain which is given by 



92 < £0 



e'n + ei 



(49) 



where 



LOfb_ 

7m 

KUJfb 



(' + ^) 

UJfb 
■ 7mK 



K7„ 



CO/(Bn 



Jfb 



7mW/6 



+ UJft (K + 7m) 



(50) 



(51) 



FIG. 3: Plot of the effective mechanical frequency of Eq. (|46p 
(a), and of the effective mechanical damping of Eq. (I47|) (b), 
versus frequency. Parameter values are LJm/27r = 10 MHz, 
7m/27r — 100 Hz, gcd = 4, = 3cJm, K = 30wm and G — 0.2 
cjm (dashed line), G = 0.3 (dotted line) and G = 0.4 ujm 
(full line), which we shall see later correspond to an optimal 
feedback cooling regime. 



the transfer function g{uj), and the one associated with 
the choice of Eq. ([45)1 is plotted in Fig. 3 for comparison 
with the corresponding curves for the back-action cooling 
of Fig. 1. Fig. 3a shows again that in the chosen parame- 
ter regime, the frequency shift, i.e., the optical spring ef- 
fect, is negligible. Cold damping is usually applied in the 
adiabatic limit when k, w/f, ^ to and in this limit one has 



(jj 



eff,cd , 



UJm and 7^f ■^'^^'^ ~ -frn+gcdGuJm/ = 7m(l+52), 

where we have defined the scaled, dimensionless feedback 
gain 92 = 9cdGujm/ Kjm In this limit the only ef- 

fect of cold damping is to increase the mechanical damp- 
ing rate without significantly changing the resonance fre- 
quency As a consequence, the system is stable whenever 
gcd > 0. This is not any longer true in the general case 
when finite values of k and lo fb are considered. By impos- 
ing that all the poles of Xejfi^) in the lower complex 



The system may become unstable for large gain because, 
for non-zero time delay, the feedback force can be out-of- 
phase with the oscillator motion and become an acceler- 
ating rather than a viscous force. 

The exact expression for the position and momentum 
variances can be obtained by integrating Eq. (fT3|) using 
the corresponding spectrum for the cold-damping case, 
given by Eq. (j42[) . Using again the approximation of 
Eq. dni), one gets 




Acd 



B. 



cd 



Gcd 



G^ 



Bed] + ^ [Bed + Gcd] > , (52) 



to 



fb 



Br 



D. 



cd 



ald'^jb 



A. 



cd 



(53) 



-^Acd + Dcd 



91 

2k 



-A 



cd 



''fb 



B. 



cd 
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where 

Acd = [lii^m + ^fb (t^m + 7mK + gcdGuOm)] , (54) 

Bed = ^mK^ {im + K + UJ fb) , (55) 



Ccd = 



UJfbK, 



Lujf, (k + Jrn) + ^ fb + 7m) 



+ 7m (t^m + + K7m)] , (56) 

Dcd = Wyj (cj^ + + K7m) + (k + 7m) ffcdCWm] 

+ ^/b ('^m + «^7m) (t^m + ^7™ + QcdGuj^) 

+ «^^m ('^m + 7mK) • (57) 

Eqs. (f52|) and (|53p show that also with cold damping 
{Sq^) 7^ (Sp^), i.e., energy equipartition does not hold in 
general. 

The optimal cooling conditions in the cold-damping 
scheme can be obtained by minimizing the sum of the 
variances of Eqs. ([5^ and ([55]) . which is non trivial in 
general. Since, however, cold-damping feedback is de- 
signed to work only within the adiabatic limit, we can 
restrict the discussion to the bad cavity limit where 
K ^ u;m,"tm- In fact, the feedback force is an addi- 
tional viscous force that is able to overdamp the me- 
chanical oscillator only when the output signal SY'^'^*{t) 
is proportional to the oscillator position dq{t), which hap- 
pens when K is much larger than the relevant frequencies 
Lu of the mechanical motion. In the good cavity limit 
(k ^ LOm, 7m) ou the contrary, the output signal SY'^'''^{t) 
is proportional to the time integral of the oscillator posi- 
tion 6q{t) and therefore the feedback force is proportional 
to the oscillator position rather than to its velocity. This 
means that in the good cavity limit the feedback loop 
has no cold-damping effect, because it increases the me- 

1 /2 

chanical frequency, uj'f^-^''^'^ ~ [w^ -I- gcdGcJm] , with- 
out appreciably modifying the mechanical damping (see 
Eqs. (gH, (gel), and (jTZl). 

We discuss the expressions of ((Jg^) and ((5p^) in the 
adiabatic limit by distinguishing two situations that de- 
pend upon the value of the feedback bandwidth ujfi,'. i) 
very large bandwidth, w/f, 3> k 3> ujm,"frm where the 
feedback is practically instantaneous, ; ii) finite band- 
width, bjfb ^ iOrn ^ 7m- In the first case one has 



1^(1+32^) 



1 + .92 
(1+52^) 



(58) 



l + <?2 



32 ^fblm 



(59) 



where we have defined the scaled dimensionless input 
power ( = 2G^/K'-fm- These results provide the general- 
ization of the results of 2^ 21 1 , where the quantum limits 
of cold-damping have been already discussed within the 
adiabatic limit. In fact, Eqs. (f58l) - (f59l) reproduce the re- 
sults of Ref. 2l| in the large-bandwidth limit of the feed- 
back except for the addition of the non-adiabatic correc- 
tion term g2jm/K for both {Sq'^') and ((5p^). The almost 



instantaneous feedback regime w/b ^ k ^ Wm, 7m is not 
convenient for cooling because of the last contribution to 
which is very large since it diverges linearly with 
ujfb- This is due to the fact that the derivative feedback 
injects a large amount of shot noise when its bandwidth 
is very large. 

In the other limit where the feedback delay time is com- 
parable to the oscillator timescales, that is, k ^ ujfb ~ 
^ 7m J one has 



(Si') 



92 



HJ)( 



1 





1+52 


+ ^ 






-1 




52 


, ^m 




'gl 


^fb 




4r/C 


1 




(■• 




" 2 







(60) 



g2lml^fb 



921 m 
UJfb 



21 



(61) 

which 
and 



These results generalize those of Refs. 
were restricted to the adiabatic limit k —^ oo, 
partially modify their conclusion that one can always 
achieve ground-state cooling in the large feedback-gain 
limit g2 ~ ^ ^ CX3 (77 = 1). This is not true in general. 
From Eq. ([50)) it is clear that, when 52 — C ^ 1; ^/b/^m: 



one obtains {Sq"^) :^ 1/2- 



-^m/4^/b+(^/52)(l+^^/^;^b), 

which implies that the ground state is approached when 
^m/^fb <C 1 and n/g2 ^ 1. However, in the same limit, 
Eq. m yields (V) ^ 1/2+ (c^2^/4c^^)(527m/c^/fc) + 
n/.g2(l + t^rn/^jb + 527m/w/b), showiug that the feed- 
back bandwidth ojfb cannot be too large, because other- 
wise ((5j)^) becomes too large. The best cooling regime 
is instead achieved for ujfb ''^ 3wm and 52 — C (i-S- 
gcd — '^G/ujm), i.e. for large hut finite feedback gain. 
This is consistent with the fact that stability imposes 
an upper bound to the feedback gain when k and ujfb 
are finite. The optimal cooling regime for cold damp- 
ing is illustrated in Fig. HI where Uef / is plotted versus 
the feedback gain g^d and the input power P, at fixed 
K = 3wm (bad-cavity condition) and ujfb = 3.5wm- We 
find a minimum value n^ff ~ 0.2, corresponding to an 
effective temperature T ~ 0.27 mK for 172 — ^ — 10^. 
Lower values of rig/ / can be obtained only if the quality 
factor LUm/jm is further increased. 



CONCLUSIONS 

We have developed a general quantum Langevin treat- 
ment of radiation-pressure ground-state cooling of a 
micromechanical oscillator, extending previous treat- 

US, 



ments 



17 



21, 28 



29| to the full parameter range 
of a stable cavity. Both cavity self cooling and cold 
damping are able to approach the ground state, and the 
comparison of the optimal cooling conditions for both 



9 




FIG. 4: (Color online) a) EflFective mean vibrational number 
Tie/ / versus the feedback gain g^d and the scaled input power 
P/Po {Po ~ 100 mW), around the optimal cooling regime for 
U3m/2TT = 10 MHz, 7m/27r = 100 Hz, m = 250 ng, a cavity of 
length L = 0.5 mm driven by a laser with wavelength 1064 
nm and bandwidth k = Suim.- The feedback bandwidth is 
ujfb = 3.5tJm. The oscillator reservoir temperature is T = 0.6 
K, corresponding to n ~ 1250. The minimum value n^// ~ 
0.2 corresponds to an effective temperature Te// — 0.27 mK. 
(b) Effective mean vibrational number n^ff versus Hi/iUm and 
uJfb/iOm at fixed input power Po = 100 mW and feedback gain 
gcd = 0.8. The other parameters are as in (a). 



schemes shows that self cooling is preferable for a good 
cavity (k < LOm)^ while cold damping is more convenient 
for a bad cavity (k > a;,„). 
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